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ON INTERNAL DAMPING OF ROTATING BEAMS 
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SUMMARY 


For a beam rotating at out a transverse axis and harmonically 
vibrating in a direction normal to its plane of rotation, it is assumed 
that an internal damping force exists inversely proportional to the 
vibration frequency and directly proportional to the time rate of 
change of the elastic load. It is shown that, with such a force, the 
internal logarithmic decrement of the beam increases asymptotically 
with the principal mode of vibration to the value ng, which it would 
have in any mode if the same beam were not rotating. The term g in 
this value is a dimensionless internal— dancing coefficient depending 
on the beam material. If the beam performs vibrations in the plane of 
rotation, then the internal logarithmic decrement in the fundamental 
mode will be either equal to, or slightly less than, ng, but will 
decrease in the second mode and will then increase asymptotically with 
the higher modes to the value ng. Thus, rotation of a beam in general 
diminishes the internal logarithmic decrements. 


INTRODUCTION 


The purpose of this investigation is to derive the implications, 
for rotating beams, of the assumption Pd = ^ (Sly")" which was 


originally made for nonrotating beams and which has been found to agree 
with experimental results for such beams. Such an investigation may be 
of particular interest in the design of helicopter blades, propellers, 
or turbine blades, which often act like beams. (See reference 1 .) 
Internal damping is especially important in flutter phenomena, where it 
increases the stability, and in vibrations of fixed— ended helicopter 
blades in the plane of rotation, where it is of the same magnitude as 
the aerodynamic load. Experimental data on the internal logarithmic 
decrements of beams rotating about a transverse axis appear to be 
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lacking, and there is therefore no conclusive empirical check at present 
on the "basic assumption used here for rotating "beams. 1 

The mechanism of internal damping in structures appears as yet to 
be not^ entirely understood. It is possible, nevertheless, to take 
internal damping mathematically into account. in vibration phenomena by 
making simple assumptions such that _ the indications will be in accord 
with experimental evidence. 

Numerous experiments (references 2 to 5 ) have been made on non- 
rotating vibrating beans, all leading to the conclusion that the 
internal logarithmic decrements of such beams are independent of the 
frequency, and therefore of - the mode, of vibration. In view of these 
results Theodorsen and Garrick (reference 6 ) have, for purposes of 
flutter Investigation, introduced a damping force in bending of the 

form p^ = ign e 2 m’y. Reissner (reference 7 ) more recently has 
suggested a damping force of the form p^ = dgg ^ (Ely")", where the 

coefficient d s is to be chosen so that the logarithmic decrements 
will be independent of the frequency . 2 

It~ will be observed that in both of the expressions for p* the 
dancing force is assumed to be in phase with the velocities, although 
the second expression does not contain an imaginary factor. It should 
be observed, moreover, that the term (Ely")" represents the unit 
elastic^ load of a rotating, as well as of a stationary, beam, while the 
term cn e 2 m’y represents exactly the unit elastic load of a stationary 
beam but only approximately^ the unit elastic load of a rotating beam. 
Although the two basic assunqptions discussed here are in this sense not 
quite equivalent for rotating beams, it will be found that in *actual 
cases they lead to virtually the same results for the logarithmic 
decrements. It is therefore necessary to treat only one of these 
assumptions in detail. 

^It may be remarked here that in a rotating beam aerodynamic as 
well as internal damping forces exist. Consequently, the effect of 
aerody nami c damping mush be taken into account in the interpretation of 
any experimental data on the logarithmic decrements of rotating beams. 

2 Previous investigators, such as Sezawa (reference 8 ), also used 
an internal damping force proportional to d/dt (Ely")", but for 
simplicity it was assumed that the proportionality factor was a constant, 
independent of the frequency of vibration. This led to results which 
were not in agreement with experiments. 

3 It can be easily shown that if the centrifugal loads had no 
effect on the modes of deflection of a rotating beam this term would 
represent the unit elastic load of a rotating beam exactly. 
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SYMBOIS 


A cross— Bectional area 

Aq cross-sectional area at root of beam. 

1 modulus of elasticity of beam material 

f _ Q cn 2 

cn “ J2~ 

f = %IL. 
en 

g dimensionless internal— damping coefficient depending on beam, 

material 

i = ^ 


I moment of inertia of structural cross section of beam 


Io 

K 

l 

m* 

n 


value of I at root of beam 


dimensionless bending-stiffness parameter 


' EI o \ 
pA oft 2 17 


length of beam 

mass per unit length of beam 

mode of vibration 


p conjilex frequency; if p = -R ± Icn (R and a> real) 

then a>/2ir is the natural frequency in. cycles per 
second, wMle 2kR/cd is the logarithmic decrement 
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<1 - p/ft 
r 
B 
t 

y = 7/1 
7 

X 

6 

I = x/l 


ft 


CD 


CO 


©n 


CO 


cn 


radial distance of- a "beam element from, axis of rotation 
radial distance of beam tip from axis of rotation 
time 

vibrational bending deflection of beam at’ any point x 
distance along beam 

internal logarithmic decrement; logarithm of_ratio of ampli- 
tude, of vibration at one time to amplitude one period 
later 

density of beam material 

dimensionless centrifugal— force parameter (r -r- t d£ 

A o l 


angular speed of rotating beam, radians per second 

natural .frequency of vibration of rotating beam, cycles 
per 2jt seconds 

natural frequency of vibration of-beam in "nth." mode when it 
is stationary (fi = 0), cycles per 2it seconds 

natural frequency of vibration of rotating beam in nth mode 
if it had no bending stiffness (E = 0), cycles per 
2 it seconds 



* = S/Sx at first; 1 = b/bi in equation (3) and thereafter 

• = b/bt 


V 


V 


¥ 
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THEORY 


Let the internal or structural damping force in a rotating, as 
veil as a stationary, vibrating beam be proportional to the elastic 
load per unit length of the beam., and let this force be in phase vith 
the vibrational velocities. Then the unit damping force p^ can. be 
expressed in the form (reference 7) 

5d = d sS |jT (W)" (1) 


where g is a dimensionless internal-damping coefficient character- 
istic of the material of the beam, and where d s is a factor which, in 
accordance with empirical data, can be chosen so that the internal 
logarithmic decrement of a stationary beam performing harmonic motion 
will be independent of the frequency of vibration. It can be shown 
that d s , thus chosen, must have the value 



( 2 ) 


(in reference 7j dg vas incorrectly expressed as d s = tu. ) 

The equation of the small free vibrations of a beam rotating with 
angular velocity about a transverse axis passing through the root 
of the beam and parallel to the vibrational displacements y of the 
beam can then be written in the following dimensionless form: 


n ( 1 y" 


a y 


+ -y ( 7T ) -Cry’)’ +^«0 


“ VI 0 n 


Ao ft* 


(3) 


where, and henceforth, * = — . The term (Ty>)' represents the unit 

oi 

centrifugal load. 

For harmonic motion, let 


y = yU)©^ 


(4) 
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Then, with the assumption (justified, a posteriori, as is shown subse- 
quently) that 1 + (pg/co) X 1 + ig, the following equation is obtained 
for y( I ) and q = p/ft : 


K(1 +■ ig)( ^r~ y" 


(Ty’ )' 



(5) 


An expression for the internal damping in any mode can be obtained 
from equation (5) in the following manner. The natural frequency of 
any- undamped mode n of - vibration, as determined by .equation (5) 
with g * 0, can be expressed, to a satisfactory approximation, by 
(reference 9) 



= -f E - f 
en on 


( 6 ) 


where f en and f cn are positive nondi mens i onal quantities repre- 
senting, respectively, the contribution of the elastic resistance (ft = 0) 
and the contribution of the centrifugal load (K = 0) to the natural 
frequency of the nth mode. 

From equation (5) it is seen that- the mathematical effect of 
internal damping here is to replace the constant K by the 
constant K(l + ig). It follows from equation (6), then, that the 
complex frequency ratio q^ for the nth mode with internal damping 
will be given by: 


<ldn 2 = + *8) 


- f , 


cn 


(T) 


Assuming, as is actually the case, that g « 1, equation (7) leads to 
the following expression for q^ n : 


'bin " 


i 


^ _ i gg; 

® 2 Won/® 


( 8 ) 


where 0) on is the natural frequency of the rotating beam in the nth 
undamped mode and is given by: 


£ 0 , 


on 


- B\l fen* 1 + f, 


cn 


(9) 
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Ob s erring that 


“en l 2 


= Kf 


en 


where co is the natural frequency of the nth mode when the beam, is 
not rotating, it is seen that equation (8) can be written as: 


a 1 °W 1 „ “Wfo 2 

^ n 2 8 <D 0n /!> 

From equation (10), one finds 



Since g « 1, it is evident that 


( 10 ) 


l + -gftjl + ig 
CO 


as was originally assumed. 

Equation (10) implies that to a first (and for practical purposes, 
sufficient) approximation, internal damping does not affect the natural 
frequency of any mode, since the imaginary part of the complex frequency 
ratio q remains unchanged by the damping. However, the damping adds 
a negative real part to the value of q, which implies the following 
value for the logarithmic decrement 6 n of any principal mode: 


5 n = *S 




en 


Kf, 


co. 


= itg 


en 


'on 


Kfen + 


-cn 


( 11 ) 


where co en 2 /co orL 2 is the ratio of the square of the natural frequency 
of the beam in any mode n when it is not rotating to the square of the 
natural frequency in that mode when it is rotating. 
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Since in general for any mode (o)^ = cn^ only 

if fl = 0), and since the ratio cD en /a> on approaches unity as the mode 
of vibration increases, it follows from equation (11) that the internal 
logarithmic decrement of a rotating "beam increases asymptotically with 
the mode to the value 5 = jtg, . which it would have in any mode if the 
"beam were not rotating. Rotation of a "beam evidently decreases the 
logarithmic decrement" in any mode of vibration of the beam, but this 
decrease will in actual cases be found appreciable mainly in the lower 
modes. As can be seen from equation (11), the decrease in the loga- 
rithmic decrements is due to the increase in the natural frequencies 
caused by the rotation of the beam. 

If the beam is rotating about a transverse axis perpendicular, 
instead of parallel, to its vibrational displacements, so that the 
vibrations are in the plane of rotation, then the basic differential 
equation remains the same as equation (3) except that the unit centrif- 
ugal load is now represented by an additional term —A/ A Q y on the left 
side of the equation. The approximate relation, equation (6), remains 
valid., although the numerical values of f cn will now be different.^ 
All of the previous reasoning leading to equations (7) to (11) thus 
remains valid and, therefore, the general equations which have been 
derived here for a rotating beam performing vibrations in a direction 
perpendicular to its plane of rotation are valid also for a beam 
vibrating in its plane of rotation. 

In the latter case, however, the value of f cn in the fundamental 
mode is usually either zero or very small. Consequently, according to 
equation (11), the internal logarithmic decrement of a beam vibrating 
in the plane of rotation will have a value of almost jtg in the 
fundamental mode, will decrease in the second mode, and will then 
increase asymptotically with the higher modes to the value rtg. 

It may be noted that all the general relations which have been 
developed here are valid regardless of the boundary conditions of the 
beam, although the numerical values of f en depend on the boundary 
conditions . 


^It can be shown that if the rotational axis passes through the 
root of the blade, then (f cn ) 2 = (fcn)i ~ wIiere the subscripts 1 
and 2 refer, respectively, to vibrations normal to, and in, the plane 
of rotation. 



NACA TN 1996 


9 


NUMERICAL EXAMPLE 


A numerical example may "be given to illustrate the nature of the 
results derived, here. Consider a uniform fixed— ended beam vibrating 
normally to its plane of rotation Then the values of f eTI and f on 
for the various inodes of vibration are given byi 

fgn 1 / 1 *' = 1.874, 4.71, 7 . 85 , . . ft(2n - l)/2 (approx.) 

n = 2 

f cn = 1 > * * *J “f 2 * 1 - 1) 

n = 1 

A typical value of K for a fixed-ended single— tubular-spar helicopter 
blade is E = 0,004. For propellers and turbine blades, the value of E 
may be higher. Figures 1 and 2 show the variation of 6 n with the 
principal mode of vibration and with the vibration frequency, respec- 
tively, for the uniform fixed-ended beam with K = 0 . 005 , K = 0 . 05 , 
and E = 1. 

If a uniform fixed-ended beam is vibrating in itB plane of 
rotation, arid if the axis of rotation passes through Its root, then the 
values of f en remain as given, but the values of f cn are now‘ 


fcn ~ 5, 14, * » *, ( 2 n + 1 ) (n 1 ) 

Figure 3 shows the variation of 5 n with the principal modes of 
vibration of such a beam. 


CONCLUSIONS 


If a beam rotating about a transverse axis is harmonically 
vibrating in a direction normal to its plane of rotation and if it is 
assumed that an internal damping force exists inversely proportional 
to the vibration frequency and directly proportional to the time rate 
of change of the elastic load, then it follows that the internal loga- 
rithmic decrement of such a beam will increase asymptotically with the 
mode of vibration and will approach the value itg which it would have 
in all modes if the beam were not rotating. The term g in this value 
is a dimensionless internal— damping coefficient depending on the beam 
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material. If the team, is Titrating in, instead or normally to, its 
plane of rotation, then the internal logarithmic decrement will he equal 
to, or will he slightly less than, the value «g in the fundamental 
mode, will decrease in the second mode, and will then increase asymp- 
totically with the higher modes to the value irg. Thus, rotation of a 
heam in general diminishes the structural logarithmic decrements. The 
greatest variation of the logarithmic decrement with the mode will in 
general occur in the lower modes . 


Polytechnic Institute of Brooklyn 

Brooklyn, N. Y. , February 1 , 19^9 
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Figure 1.- Variation of structural logarithmic decrements with principal mode 
of vibration for a uniform fixed -ended beam rotating at angular speed n 
about a transverse axis and vibrating normally to the plane of rotation. 
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Figure 2.- Variation of structural logarithmic decrements with vibration 
frequency for a uniform fixed -ended beam rotating at angular speed n 
about a transverse axis and vibrating normally to the plane of rotation. 
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Figure 3.- Variation of internal logarithmic decrements with principal mode 
of vibration for a uniform fixed-ended beam rotating at angular speed n 
about a transverse axis, through its root and vibrating in the plane of 
rotation. 
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